Semiclassical Mechanism for the Quantum Decay in Open Chaotic Systems 
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We address the decay in open chaotic quantum systems and calculate semiclassical corrections to 
the classical exponential decay. We confirm random matrix predictions and, going beyond, calculate 
Ehrenfest time eflFects. To support our results we perform extensive numerical simulations. Within 
our approach we show that certain (previously unnoticed) pairs of interfering, correlated classical 
trajectories are of vital importance. They also provide the dynamical mechanism for related phe- 
nomena such as photo-ionization and -dissociation, for which we compute cross section correlations. 
Moreover, these orbits allow us to establish a semiclassical version of the continuity equation. 
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Besides their relevance to many areas of physics, open 
quantum systems play an outstanding role in gaining an 
improved understanding of the relation between classical 
and quantum physics [Ij]. For a closed quantum system 
the spatially integrated probability density 



P{t) 



dvi,{v,t)r{v,t) 
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of a wave function '0(r, t) in the volume V is constant, 
i.e. p{t) = l. This fact is naturally retained when taking 
the classical limit in a semiclassical evaluation of Eq. ([1]), 
reflecting particle conservation in the quantum and clas- 
sical limit. However, when opening up the system, p{t), 
then representing the quantum survival probability, ex- 
hibits deviations from its classical counterpart Pc\{t)] in 
other words, certain quantum properties of the closed 
system can be unveiled upon opening it. 

For an open quantum system with a classically chaotic 
counterpart, the classical survival probability is asymp- 
totically Pc\{t) = exp(— t/rd), with classical dwell time 
Trf. This has been observed in various disciplines, either 
directly, as in atom billiards 0, or indirectly in the 
spectral regime of Ericson fluctuations in electron Q or 
microwave cavities, and in atomic photo-ionization Q. 

However, it was found numerically l7| and confirmed 
with supersymmetry techniques [^^that the differ- 
ence between p{t) and pc\{t) becomes significant at times 
close to the quantum relaxation time t* = ^JtsXh. In 
the semiclassical limit t* is shorter than the Heisenberg 
time Ih — l-Khj IS. (with A the mean level spacing). It 
was shown in Refs. [1, that pit) is a universal func- 
tion depending only on and tu in the random matrix 
theory (RMT) limit Though the leading quantum 

deviations from Pq\{^) were reproduced semiclassically for 
quantum graphs ll|, a general understanding of its dy- 
namical origin is still lacking. 

In this Letter we present a semiclassical calculation of 
pit) for t < t* for general, classically chaotic systems. 
It reveals the mechanism underlying the appearance of 
quantum corrections upon opening the system. In our 



calculation we go beyond the so-called diagonal approx- 
imation and evaluate contributions from correlated tra- 
jectory pairs IJ]. This technique has been extended and 
applied to calculate various spectral [3, [3] and scat- 



tering 15, 3, 17, 3, 3 13, 21 1 properties of quantum 
chaotic systems. We find however that for calculating 
p{t) a new class of correlated trajectory pairs, 'one-leg- 
loops', has to be considered along with the previously 
known loop diagrams. They prove particularly crucial 
for ensuring unitarity in problems involving semiclassi- 
cal propagation along open trajectories inside a system 
and, moreover, allow one to semiclassically recover the 
continuity equation. 

We present the dominant quantum corrections to Pci(t) 
for systems with and without time reversal symmetry. 
Going beyond RMT, we calculate Ehrenfest time effects 
on p{t), compare with quantum simulations of billiard 
dynamics, and extend our approach to photo-ionization 
and -dissociation cross sections. 

Semiclassical approach - We consider p{t), Eq. ([1]), for 
a two-dimensional system of area A and express 4^(v, t) = 
J dr' K{r,r' ■,t)tpQ{r') through the initial wave function 
-00 (r') and the time-dependent propagator K{r, r'; t) that 
we approximate semiclassically by 22 1 



Ksc (r,r';i) 
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7(r'^r,t) 



Here = Sj{r,r';t) is the classical action along the 
path 7 connecting r' and r in time t, and — 

1 1 /2 

\det{d'^ / drdr')\ exp(— i7r/i;y/2) with Maslov index p^. 

The semiclassical survival probability, psc (t) , obtained 
by expressing the time evolution of ip{r, t) and ip*{r, t) in 
Eq. ID) through K^c^ Eq. is given by three spatial in- 
tegrals over a double sum over trajectories 7,7' starting 
at initial points r' and r", weighted by V'o(r') and i/'o (>""); 
and ending at the same point r inside A. For simplicity of 
presentation we here assume ipo to be spatially localized 
(e.g. a Gaussian wave packet); while generalizations are 
given below. Introducing ro = (r'+r")/2 and q= (r'— r"). 
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we replace the original paths 7,7', by nearby trajecto- 
ries 7 and 7' connecting Tq and r in time t. Then, upon 
expanding the action S^{r,r';t) ~ 5^(r,ro;i) — qPo/2 
(with Pq the initial momentum of path 7) and S^i anal- 
ogously, we obtain 

Psc(i) = j (^^^^odciiJoiro + |) (i"o - |) 

X J2 i^^^;,e(*/'')[^--V-(pS+pS')q/2]. (3) 

7.7'(ro— »r,t) 

The double sums in Eq. ^ contain rapidly oscillating 
phases {S-y~Sryi)/h which are assumed to vanish unless 7 
and 7' are correlated. The main, diagonal, contribution 
to psc arises from pairs 7 = 7' that, upon employing a sum 
rule [ill, yield the classical decay pci(i) = {e-*/^") (H]. 
Here {F) = {2iTh)-'^ j drodpo-F(ro, Po)pw(ro, Po), where 
PH'(ro,Po) = / dqV'o(roW2)^o*(ro-q/2)e-(*/'"')'i-P« is the 
Wigner transform of the initial state and Td = il{E)/2wp, 
with U{E) = J drdp6{E — H{r,p)) and w the size of the 
opening. For chaotic billiards Tcj(p) =m7r A/ wp. For ini- 
tial states with small energy dispersion pci{t)—e^*^'^''^P"\ 
For systems with time reversal symmetry, leading- 
order quantum corrections to Pci{t) arise from off- 
diagonal contributions to the double sum in Eq. ([3]) , given 
by pairs of correlated orbits depicted as full and dashed 
related semiclassical treatments 
The two orbits are expo- 
nentially close to each other along the two open 'legs' 
and along the loop [l5| . but deviate in the intermedi- 
ate encounter region (box in Fig.[Ha)). Its length is 
icnc = A~"'^ln(c^/|su|) [l3|, where A is the Lyapunov expo- 
nent, c is a classical constant, and s and u are the stable 
and unstable coordinates in a Poincare surface of section 
(PSS) in the encounter region. Such 'two-leg-loops' (211) 
are based on orbit pairs with S-y — Sji = su and a density 
W2ii(s,M,i) = [t- 2ic„c(s,M)]V[217(£;)tc„c(s,w)][23l. In- 
voking the sum rule, the double sum in Eq. ([3]) is replaced 
hy J du J dse(-*+*-<=)/^'^w;2ii(s,M,i)e('/'*)"". Here e*"-/^'' 
accounts for the fact that if the first encounter stretch is 
inside A the second must also be inside A. This gives the 
211 contribution (Fig.IlJa)) to p{t): 



P2\\{t) 
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The linear term in Eq. ^ violates unitarity, since it does 
not vanish upon closing the system, i.e. as — > 00. 
This is cured by considering a new type of diagrams. 
These orbit pairs, to which we refer as 'one-leg-loops' 
(111), are characterized by an initial or final point in- 
side the encounter region (Fig.[ljb,c)). They are relevant 
for open orbits starting or ending inside A and hence 
have not arose in conductance treatments based on lead- 
connecting paths, since at an opening the exit of one 
encounter stretch implies the exit of the other one 24 1. 
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FIG. 1: (Color online) Pairs of correlated classical trajectories 
7 (full line) and 7' (dashed) generating the leading quantum 
corrections to the classical decay probability. While in panel 
(a) the encounter region (box) connects a loop with two legs, 
the paths begin or end inside the encounter region ('one-leg- 
loops') in (b) with and in (c) without a self-crossing in con- 
figuration space. The zoom into the encounter region in (b) 
depicts the position of the Poincare surface of section used. 



For their evaluation consider the time t' between the 
initial or final point of the trajectory and the PSS, de- 
fined in the zoom into Fig.Ilfb). Then teiic{t',u) = 
— t' + A^^ln(c/|M|) and Sj — S^' su for any position 
of the PSS. The density of encounters is wii\{s,u,t) = 
2 jA"M„(c/|a|) _ 2t,^^{t',u)]/[n{E%nc{t',u)], where 
the prefactor 2 accounts for the two cases of beginning 
or ending in an encounter region. We evaluate this con- 
tribution by modifying Pci(i) by e*"'"/'^'' as before and 
integrating over s, u and t' . To this end we substi- 
tute [H t" = t' + ln(c/|u|), a = c/u and x = su/c^, 
with integration domains — l<x<l, I<cr< e'^* and 
0<t" <\-^ln{l/\x\). Note that the limits for t' include 
the case when the paths do not have a self-crossing in 
configuration space (Fig.[TJc)). The integration yields 



Pill (t) 



2±e 

tH 



-t/Ta 



(5) 



It precisely cancels the linear term in p2i\, Eq. i.e. 
P2ii(<) + piii(t) = e~*/^''t^/(2Tdt/f), recovering unitarity! 

The next-order quantum corrections are obtained by 
calculating 111 and 211 contributions of diagrams such 
as discussed in [l^. Together with Eqs. (|4l5p . this yields 
for systems with time reversal symmetry 



Psc(t) 



5t^ 



2Tdt 



H 



H 



(6) 



for t < t* . The term quadratic in t represents the weak 
localization-type enhancement of the quantum survival 
probability. The quadratic and the quartic terms, which 
agree with RMT [8| , dominate for the time range consid- 
ered. The cubic term in Eq. ([S]), whose functional form 
was anticipated in scales differently with and tn. 
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For systems without time reversal symmetry the calcu- 
lation of the relevant one- and two-leg-loops gives, again 
in accordance with RMT m , 



Pscit) 



1 



(7) 



We finally note that our restriction to localized initial 
states can be lifted and the results generalized to arbi- 
trary initial states 25|. This is because the trajectory 



pairs 7, 7' survive the integration in Eq. ([3|) only if their 
starting points are close to each other in phase space, 
rendering the above analysis valid. 

Continuity equation.- It is instructive to reformulate 
the decay problem in terms of paths crossing the open- 
ing. To this end we consider the integral version of the 
continuity equation, 9p(r, <)/9i+V-j(r, i) = 0, namely 



d f 



(8) 



where S is the cross section of the opening with a nor- 
mal vector fix- In Eq. ([8]), the current density j(r, i) — 
(l/m)Re[(?i/i)'(/'*(r, t)V'(/'(r, i)] can be semiclassically ex- 
pressed through Eq. ([2]) in terms of orbit pairs connect- 
ing points inside A with the opening. In the diagonal ap- 



proximation we obtain /ojc 



■fixdx 



-t/ra 



/Td, consis- 



tent with Pciit). Loop contributions are calculated anal- 
ogously to those of psc from Eq. giving 



(j2n+jiii) • nxdx = e ^^^^ 



(9) 



Time integration of Eq. ([5]) leads to P2ii{t) + Pm{t) ~ 
g-t/rrf^2/^2Td%), consistent with Eq. ®. The 111 con- 
tributions enter into Eq. ^ with half the weight, since 
Ills with a short leg (encounter box) at the opening must 
be excluded. These 'missing' paths assure the correct 
form of quantum deviations from Pciit). 

Higher 211 and 111 corrections to j lead to Eqs. (|6|7p . We 
conclude that both, 211 and 111 contributions to j are es- 
sential to achieve a unitary flow and thereby to establish 
a semiclassical version of the continuity equation. 

Ehrenfest time effects. - The Ehrenfest time te (26j sep- 
arates the evolution of wave packets following essentially 
the classical dynamics from longer time scales dominated 
by wave interference. While r^; -effects have been mainly 
considered for stationary processes involving time inte- 
gration [3, [13, H [Slil [ill, signatures of te should 
appear most directly in the time domain 3, 2^, i.e. for 
p{t). Here we semiclassically compute the r£;-dependence 
of the weak-localization correction to p{t) in Eq. ©. 
To this end we distinguish between t|; = A^^ln(£/A_B), 
where C is the typical system size and Xb the de Broghe 
wavelength, and = X~^\ii{'w'^ / (CXb)) , related to the 
width w of the opening [l^. As before we consider that 
the densities W2ii,m{s,u,t) contain the Heaviside func- 
tion 9{t — 2tcnc) (negligible for t'^'^ <^ Td) assuring that 




t/TH 



FIG. 2: (Color online) Averaged ratio R(t) between numerical 
quantum mechanical decay pqm (*) and classical decay Pci™(i) 
(red symbols, the bars correspond to the variance after av- 
eraging, see text) compared with corresponding semiclassical 
predictions based on the quadratic term in Eq. ((6} (dashed 
line) and Eq. JTOj (full line). Upper inset: Pqm(t) (red full 
line) and p^J™(f) (dashed) for the wave packet shown in the 
lower inset. Lower inset: Desymmetrized diamond billiard, 
defined as the fundamental domain of the area confined by 
four intersecting disks of radius R centered at the vertices 
of a square of length 2L (L — 100, R — 131) with opening 
w — 16. The initial Gaussian wave packet shown is of size 
a — 3 and Ab = 3. The arrow marks the momentum direction. 



the time required to form a 111 or 211 is larger than 2tonc- 
Our calculation gives (for t'^'^X ^ 1 with 2t'^—t'^ + t^) 



P2\i{t) + Pn\it) = e 



2Tdt 



2r2 



H 



(10) 



Numerical simulation.- The leading-order quantum 
corrections in Eq. ^ and ([7]) were confirmed by nu- 
merical simulations for graphs 11]. Here we compare 
our semiclassical predictions with quantum calculations 
of p{t) based on the numerical propagation of Gaussian 
wave packets inside a billiard, a setup much closer to ex- 
periment. We chose the desymmetrized diamond billiard 
(inset Fig. [2]) [13] that is classically chaotic (A^^ ~ 3t/, 
with Tf the mean free flight time). Its opening w corre- 
sponds to iV = 10 open channels and to ~ 15''"/ • For 
the simulations, we reach t^ ~ 10. 6r^ implying t* = 3.3Td, 
T° ~0.17Trf and r|;~0.55r<j (with C^^/l). 

In the upper inset of Fig. [2] we compare the decay, 
Pq™(t) (red full line), for a representative wave packet 
simulation with the corresponding classical, pfi^{t) 
(dashed line), obtained from an ensemble of trajectories 
with the same phase space distribution as the Wigner 
function of the initial quantum state. pf\^{t) merges 
into the exponential decay exp(— i/rd), and Pq™(i) co- 
incides with pf{^{t) up to scales of t*. For a detailed 
analysis of the quantum deviations we consider the ratio 
Rit) = [ptmit)-pT'it)]/pT'it)- The red symbols in Fig.H 
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represent an average of R{t) over 27 different opening 
positions and initial momentum directions. The dashed 
and full curve depict the semiclassical results based on 
the quadratic term in Eq. ^ (dominant for the t/r^- 
range displayed) and on Eq. pH)) . The overall agreement 
of the numerical data with the full curve indicates te- 
signatures. We note, however, that we cannot rule out 
other non- universal effects (e.g. due to scars short 
orbits, diffraction or fluctuations of the effective 01) 
that may also yield time shifts. Furthermore the indi- 
vidual numerical traces R{t) exhibit strong fluctuations 
(reflected in a large variance in Fig. [2]). A numerical con- 
firmation of the log(l/fi,)-dependence of te seems to date 
impossible for billiards. 

Photo-ionization and -dissociation cross sections. - Ke- 
lated to the decay problem are photo absorption pro- 
cesses where a molecule 31| (or correspondingly an atom) 
is excited into a classically chaotic, subsequently decay- 
ing resonant state. In dipole approximation, the photo- 
dissociation cross section of the molecule excited from the 
ground state \g) is cr(e) = ImTr{AG(e)J-, where G{e) is the 
retarded molecule Green function, A = [e / {cheQ)]\(t)) 
and \(j)) = D\g), with D = d-e the projection of the dipole 
moment on the light polarization e. The two-point cor- 
relation function of a{e) is defined as 

C(c^) = {a{e + LuA/2)a{e - ujA/2))J{a{e)f^ - 1 , (11) 

where (cr(e))e w Tr{2Trh)~^ f drdpAw{r,p)6{e- H{r,p)) 
semiclassically, with Wigner transform of A. Previ- 



ous semiclassical treatments of C{uj) S^, 33| were limited 
to the diagonal approximation. To compute off-diagonal 
(loop) terms we consider Z{t)= j'^^duje'^^'^'^'^ C{lo) with 
T = t/tH- Semiclassically, Zsc is again given by a dou- 
ble sum over orbits with different initial and final points. 
Due to rapidly oscillating phases from the action differ- 
ences, only two possible configurations of those points 
contribute : (i) orbits in a sum similar to Eq. Q 
leading to a contribution as for psc{t); (ii) trajectories 
in the vicinity of a periodic orbit. Expanding around 
it, as in [s^l, leads to the spectral form factor K°^™{t) 
of an open system. From (i) and (ii) we have Zsc{t) = 
K°^°'^{t) + 2psc{t) for the time reversal case. Up to sec- 
ond order in r > we find if °p<'"(t) =e-^^(2T- 2t2) and 



Psc(t) - e-^-(l + 7VrV2) (Eq. ©), Thereby Z,, 
g-iVTj2_|_2r-|-(A^— 2)t^], confirming a conjecture of 
Its inverse Fourier transform yields the two-point corre- 
lation (with r = 2'KLOTd/tH) 



c.c(r): 



1 



ivi+r2 



1- 



1 i-r^ N-2 i-sr^ 



N i+r2 iv2 [i+v^Y 



(12) 



The first two diagonal terms agree with [34I; the third 
term represents the leading quantum correction. 

To conclude, we presented a general semiclassical ap- 
proach to the problems of quantum decay and photo 
cross-section statistics in open chaotic quantum systems. 
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